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The geometrically different plan forms of near wall plume structure in turbulent natural 
convection, visualised by driving the convection using concentration differences across 
a membrane, are shown to have a common multifractal spectrum of singularities for 
Rayleigh numbers in the range 10^° — 10^^ at Schmidt number of 602. The scaling is 
seen for a length scale range of 2^ and is independent of the Rayleigh number, the flux, 
the strength and nature of the large scale flow, and the aspect ratio. Similar scaling is 
observed for the plume structures obtained in the presence of a weak flow across the 
membrane. This common non trivial spatial scaling is proposed to be due to the same 
underlying generating process of the near wall plume structures. 



1. Introduction 

Turbulent Rayleigh - Benard convection and other related natural convection flows at 
high Rayl eigh numbers have many unresolved issues; these are extensively discussed in the 



review of Sigeia ( 199 4|). As most high Ra studi es are conducted in cryogenic conditions 
ijNiemela. Skrbek. Sreenivasan fc DonellvlEinnI) . which prevent visualisation of the near 
wall structures, little is known about the nature of near wall coherent structures in 
high Rayleigh number turbulent free convection. The non-dimensional parameters that 
characterise natural convection flows are the Rayleigh number {Ra — g^H^ /i/k, the 
ratio of buoyancy effects to dissipative effects); Prandtl number {Pr = v/k , a, fluid 
property); Nusselt number ( Nu = ^at /h ' ^ non-dimensional flux) and aspect ratio 
(AR = L/H , a geometric parameter of the fluid layer) where, g = the acceleration due 
to gravity, ^ = nondimensional density difference, Q = the kinematic heat flux, AT 
= the temperature difference between the walls, H = fluid layer height, L = fluid layer 
width, v— kinematic viscosity and k = thermal diffusivity. 

Most of the previous studies on near wall coherent structures give only a qualitative pic- 
ture of closed polygonal forms of line plumes at Ra ~ 10^, that change to randomlv mov- 
ing a nd merging line plume patterns at higher Ra of 10^ — 10^ fepangenber g fc Rowlandl 
Il96lt FKmai fc Asacda 1984; Thecrthan fc Arakeri 2000; Kerr fc Hcrrinai l200nl) . Line 
plumes are buoyant fluid rising in the form of sheets from li nes on the horizontal heated 
surface. Visualisations at the highest Ra so far (~ 10^) bv IZocchi. Moses fc LibchabeJ 
l)l990|) sho w lin e pl umes being swept along the direction of shear caused by the large 
scale flow. iTheerthan fc Arakeril l)l998(l showed that the randomly spaced and oriented 
plumes near the wall can be modeled as a regular array of laminar line plumes, each 
plume fed by laminar natural convection boundary layers on either side. IPuthenvcettiil 
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Figure 1. (a) The experimental set up (b) The plan form view of the near wall plume structure 
as seen at the intersection of a horizontal laser sheet, (c) The side view in a vertical plane 
showing rising line plumes. Planform in (b) is for the TF-type convection showing a plume free 
area where the large scale flow impinges on the membrane. 



I|2004|) extended the analysis to high Ra- high Pr case, and found that the probabihty 
distribution of plume spacings have a common log-normal form, independent of the pa- 
rameter values. None of these studies address the spatial scaling of the plume structure 
pattern. 

In this paper, we report the spatial scaling of the near wall plan-form plume struc- 
tures in high Rayleigh number turbulent free-convection, using multifractal formalism. 
Briefly, multifractal formalism describes the statistical properties of singular measures 
in terms of the singularity spectrum /(a), corresponding to the singularity strength a. 
f{a) can be regarded as the fractal dimensions of subsets with cor responding singular - 
ity strength a. The details of multifractal formalism can be seen in iMandelbrotI lll98<i): 
iMeneveau fc Sreenivasa ^ lll99j) : lHaise7' Jensen. Kadanoff. Procaccia fc ShraimanI l)l986|) 
and the references cited therein. Even though studies abound on applying this formalism 
for characterising spatial structures in various fields, we are not aware of such an analysis 
of near wall coherent structures of turbulent convection. 



2. Experimental setup and image preprocessing 

We visualise the near wall plume structure at high Ra by driving the convection us- 
ing concentration differences of NaCl across a membrane. Concentration in the present 
experiments is equivalent to temperature in Rayleigh - Benard convection. High Ra 
{10^^ — 10^^) are achieved at Schmidt number {Sc ^ v/D, equivalent to Pr) of 602 due 
to the low molecular diffusivity of NaCl. The set-up consists of two glass compartments 
of square cross section, arranged one on top of the other with a fine membrane fixed 
horizontally in between them; the schematic is shown in Figure |21 The membranes used 
are Pall Gelmann™NX29325 membrane disc filters (PG henceforth) with a random pore 
structure having a mean pore size of 0.45/x and Swedish Nylobolt 140S screen printing 
membrane (140s henceforth) with a regular square pore of size 35/i x 35//. 

The bottom tank is filled with distilled water tagged with a small amount of Sodium 
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Fluoresceine (absorption spectra peak at 488nm and emission spectra peak at 518nm) 
and then the top tank fiUed with brine to initiate the experiment. The convection is 
unsteady, but quasi-steady approximation can be used as the time scale of change of 

A 

the density difference (-^/-^f-) is 10 times the time scale of one large scale circulation 

{H/W^), where - ro/?Og^^^'^ l(Deardorfjll97(il) is an estimate of the large scale flow 
strength, H is the top tank fluid layer height, Q denotes the flux of NaCl and /3 is 
the coefficient of salinity. A horizontal laser sheet, expanded and collimated from a 5W 
Spectra Physics Stabilite™2017 Ar- Ion laser at 488nm is passed just above (<1 mm) 
the membrane. The dye in the bottom solution while convecting upward fluoresces on 
incidence of the laser beam to make the plume structure visible. A schematic of the 
visualisation of the near wall phenomena is shown in Figure |2Ib) and (c). A visible long 
pass filter glass. Coherent optics OG-515, is used to block any scattered laser light and 
allow the emitted fluorescence to pass through. The images are captured on a digital 
handycam Sony DCR PC9E. Experiments are conducted in 23 cm high tanks, with one 
tank having 15 cmxl5 cm (AR = 0.65) cross section and another with 10 cmxlO cm 
(AR — 0.435) cross section. Starting top tank NaCl concentrations of lOg/1, 7g/l and 
3g/l are used to study the plume structure under different Ra. 

The Laser Induced Fluorescence (LIF) images are RGB 24bit color with 640 x 480 
pixels resolution. The multifractal analysis is conducted on binary images obtained from 
these RGB images. In using binary images, we are neglecting the intensity variation of 
the fluorescence (proportional to the concentration of the dye) within the plume line 
thickness. The analysis is hence valid only for the geometrical aspects of the planform. A 
similar telegraph approximation of the temperature time trace measured at the middle 
of the fluid layer has been studied bv iBershadskii. Niemela. Praskovskv fc SreenivasanI 
l)2004|) to show the clusterisation of plumes. The RGB image is cropped to remove the 
tank walls, converted to gray scale. Radon transformed to remove the lines (formed due to 
imperfections in the optics and the test section surface) and then re-sampled to increase 
the resolution. The non-uniform background illumination due to the attenuation of the 
laser sheet is subtracted to obtain the plume lines over a uniform dark background. 
Contrast and intensity are enhanced to make the plume lines clearer before converting 
to a black and white binary image using a threshold. The effect of thresholding on the 
multifractal exponents is considered in Appendix The binary image retains almost all 
the features of the raw image. 



3. The plume structure planforms 

Figm-e El shows the binary images of the various planform plume structures obtained 
in the experiments. The parameter values are shown below each image. AC(g/l) is the 
effective driving concentration difference on the side of the membrane where the structure 
is visualised; Ra is based on this AC. The white lines in the figures, of thickness ~ Imm, 
are the bases of the sheet plumes originating from the membrane surface. Two types 
of convection are observed depending on the pore size of the membrane. In one case, 
the transport of salt through the membrane is by diffusion (D-type henceforth). In the 
second case, the transport is due to a flow, albeit weak, across the membrane (TF type 
henceforth). In both cases, the line plumes emanate from the thin unstable boundary 
layers above the membrane. As the focus of this paper is on the spatial scaling of these 
structures, we give only a brief description of the phenomena in each image of Figure El 
the details are discussed in .Piithenveettil (2004') . 
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(a) PG membrane in 
15 cmxl5 cm tank, 
Ra = 2.04x10", AC 
= 3.2 g/1, Q = 0.1 
mg/cm^/min, = 
0.31 cm/s, 7?e=797, Im- 
age area = 14.65cm 
square, 714^pix 



(b) Structure after 
30seconds from (a) 
Ra = 2.03x10", AC 
= 3.19 g/1, Q = 0.1 
mg/cm^/min, Wt = 
0.31 cm/s, Re = 797, 
Image area = 13.74 cm 
square, 728^ pix 



(c) PG membrane in 
10 cmxlO cm tank, 
Ra = 2.03x10", AC 
= 3.198 g/l, Q = 
0.11 mg/cm^/min, W» 
= 0.31cm/s, Re = 797, 
Image area = 9.157 cm 
square, 728^ pix 





(d) PG membrane in 
15 cmxl5 cm tank, Ra 
= 5.56x10"', AC = 
0.875 g/1, Q = 0.02 
mg/cm /min, Wr. = 
0.18cm/s, Re = 463.6, 
Image area — 12.83cm 
square, 728^ pbc 



(e) 140s membrane in 
15 cmxl5 cm tank, 
Ra = 5x10", AC = 
7.87 g/1, Q = 0.39 
mg/cm /min, W* = 
0.47 cm/s. Re = 1210.5, 
Vi = 0.002 cm/s. Image 
area = 8.61 cm square, 
528^ pix 



(f) 140s membrane in 
15 cmxl5 cm tank, 
Ra = 4.49x10", AC 
= 7.06 g/1, Q = 0.41 
mg/cm^/min, W* — 
0.48 cm/s. Re = 1236.3, 
Vi = 0.0023 cm/s. Im- 
age area = 8.32 cm 
square, 520^ pix 




(g) 140s membrane in 
15 cm X 15 cm tank, 
Ra = 5.36x10", AC 
= 8.42g/l, Q = 0.41 
mg/cm^/min, W* = 
0.48 cm/s, Re = 1236.3, 
Vi = 0.0016 cm/s, Im- 



(h) 140s membrane in 15 
cmxl5 cm tank, Ra = 
1.65x10", AC = 2.6 g/1, 
Q = 0.07 mg/cm^/min, 
= 0.27cm/s, Re = 
695.4, Image area= 14.69 
cm square, 714^ pix 
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3.1. Planforms in D-Type convection 

For the lower pore sized (PG 0.45/i) membrane, the transport across the partition is 
diffusion dominated, while the transport above and below the membrane is similar to 
Rayleigh - Benard convection at high Ra. Figures 2(a) to 2(d) show the structures ob- 
tained in the D-type convection. Similar plumes cover the cross section on the lower 
side of the membrane. The first two images are from the larger AR =0.65 which had 
multiple large scale flow cells, the signatures of which are seen as circular patches with 
aligned line plumes oriented radially around the patches, the plume free circular patch 
being the area where the large scale circulation impinges on the membrane. The diffusion 
layer thickness 5 = DAC/Q is about 0.26 mm for the image of Figure 4(a). The large 
scale flow cells shift their position at a larger time scale (~ 20 min) than one large scale 
circulation time (~ 1 min), which is greater than the time scale of few seconds for the 
merging of plumes. Figure 2(c) shows the planform of plume structure under the same 
conditions as in Figures 2(a) and 2(b) but in the smaller (10 cm x 10 cm, AR=0.435) 



cross sectional area tank. In this case, there was only a single large scale circulation 
which was rotating clockwise in the x — z plane and aligning the near wall plumes in the 
x-direction (See Figure 2(d) for the co-ordinate directions, z is normal to the membrane). 
The plume structure at lower Ra in the larger AR =0.65 tank, shown in Figure |2(d)[ 
had two counter rotating circulation cells ( anticlockwise on the left and clockwise on the 
right) rotating in the x — z plane which created a near wall mean shear directed toward 
the center along the x direction. In all these images, the flux scaled approximately as 
AC^/'^, equivalent to Nu Ra}/^ in turbulent Rayleigh - Benard convection. 



3.2. Planforms in TF-Type convection 

Experiments with the coarser (140s, 35/i) membrane showed that only half the area 
above the membrane was covered by plumes, with the other half having plumes below the 
membrane. This structure is due to a weak through flow across the membrane. The lower 
fluid dynamic resistance of the coarser membrane allows this through flow. The region 
having plumes on top of the membrane corresponds to upward through flow and vice 
versa ( Figure El^b)). The through flow velocities Vj are about 10 times smaller than the 
near wall velocity scales in turbulent free convection given by Wo ~ (aBOD)^^^l TownsenS 
I1959I) . Figure 2(e) to 2(g) show the plan form view of the plume occupied region in the TF 



type convection at AR =0.65. Figure 2(e) shows the corner region of this structure when 
there was a near wall mean shear toward the bottom left corner due to a diagonally 
oriented large scale circulation rotating in the clockwise direction. Figure 2(f) shows 
a similar corner view of the structure when there is a near wall mean shear along the 
diagonal toward the top right corner, created by an anticlockwise circulation. The central 
zoomed view, where the mean shear effects are predominant, in a structure similar to 



that in Figures 2(e) and 2(f) is shown in Figure 2(g) In all these cases, the flux scaled 
as AC^ due to the presence of a flow across th e mem brane. The phenomenology behind 
this flux scaling is described in IPuthenveettill l|2004j) . At lower driving potentials, the 



TF-type convection in the 140s membrane experiments was seen to change to the D-type 
with Q ~ AC^/'^ scaling. Figure 2(h) shows one such type of plume structure. 

In both the types of convection, the plume structures are continuously evolving spatio- 
temporal patterns. But, as the time scale of the merging of plumes is much smaller than 
the time scale of change of Rayleigh number, the planforms could be expected to exhibit 
statistically stationary characteristics for a given set of parameters. In addition, the scale 
free, non uniform and dendritic appearance, along with the the common probability 
distribution of the spacings, motivated us to analyse the structures using multifractal 
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4. Multifractal analysis 

4.1. Methodology and Scaling range 

We use the standard box counting methodology to estimate the multifractal exponents. 
The analysis is done on square binary images (L^ pixels). The measure Pij where, i, j are 
the box indices, is the probability of occurrence of the plume in a box, calculated as the 
ratio of white pixels in each box to the total number of white pixels in the image. For each 
moment q, the partition function Zq = ^ ■ ■ Pf- is calculated for various box sizes {lb in 
pixels). The slope of InZ^ Vs In (^) gives the Cumulant generating function T{q) defined 
through ■ P^- ^ i^y^'^'' where, T{q) — {q—l)Dq with Dq referring to g*'* order Renyi 
dimension. The number of boxes where Pij ^ (^) has a singularity strength between 

a and a + da is given by N{a) ~ (^) f{a) is the fractal dimension of the subset 

- picked out by each value of the moment q - with singularity strengt h a. We calculate 
the Holder exponent a and f{a) using the direct method due to IChabra fc JensenI 

l)l989 f). The expressions for /(a) and a are f{q) = lim/j^^o — TS^lT/iy^"^' ^(l) = 
limij,_,o ^ ' iiffi J/L) ^' ^ , where /i^j = ^ are the normalized measures. For each value 

of the moment q, the slope of the of the plots of the numerator vs the denominator in 
the expressions for f(q) and a(q), in the range of values where the plots are linear, give 
/(g) and a{q). In our work, the range of q values is limited to — 2 ^ g ^ 5 in order to 
have a reasonable range of scaling regime. 

The estimation of /(a) and a for the image in Figure 2(e) is shown in Figure 13 The 
image size is 8.61 cm square sampled at 528^ pixels (L= 528). The calculation is done for 
box sizes k = 528, 264, 176, 132, 88, 66, 48, 44, 33, 24 and 22 pixels. Figure |3(a)1 shows 
the plots of J fiij In /i^ j versus ln{li,/L) for a few typical values of g = —2, 0.5 and 5. 
The slope of the linear fit in the range 66 pix ^ Zf, ^ 528 pix, (i.e., from 1.076 to 8.61 cm) 
is used to estimate /(a(g)) for all q. This slope is vahd for a larger length scale range of 
22 pix ^ Zb ^ 528 pix, (i.e. from 0.359 to 8.61 cm ) for g =-0.5 and 5.0, while for g = -2, 
the linearity holds only in the fit range. The same linearity ranges for various q values 
are used for estimating a(g) from the plot of J2i j Mij Pij vs ln(/b/L) (Figure 



3(b) I 



and T(g) from the plots of In^- ^ P^^j vs \n{lb/L) (Figure 3(c) I. Figure 3(a) to 3(c) show 
that the deviation from the linearity is pronounced for g = — 2 below a specific box 
size of around 1cm (hi(lh/L) = —2), while it is almost negligible for positive moments. 
This b eha vior is also seen in other studies l ike Chabra. Mcncvcau. Jensen & SrccnivasaJ 
l)l989j) and lMeneveau fc SreenivasanI f 199l(l . As negative moments amplify low measures, 
and the effect of noise in the low measure regions at smaller box sizes is substantial, we 
expect this deviation to be due to the amplification of the errors. For the same reason, 
the negative moment slopes were also the most affected when the binary image threshold 
was changed (See inset of Figure Ej). 

The analysis is carried out for all the images in Figure [21 The major common obser- 
vations from the plots of J^i j l^i,j Mi. j J2i j Mi j Pi.j vs ln(4 /L) for these images 
are as follows, (a) For moments g ^ — 1, the linearity of the above quantities is seen 
to hold till about 0.5 cm. As positive moments pick out areas of higher measure, these 
moments represent the main plume structure. Thus, the main plume structure shows a 
multifractal scaling in the range of tank cross section to 0.5 cm i.e., a range of 2^. (b) 
When negative moments are less than -1, the plots of these quantities show a deviation 
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(a) Range for f{a) 



(b) Range for a 



(c) Range for r(g) 



Figure 3. The range of multifractal scaling for moments q — -2, 0.5 and 5. The lines show the 

linear fit used for calculating the slopes. 




Figure 4. The probability distribution function of the measure Pij in the standardised form at 
three box sizes. The inset shows the normal probability plot of the measure at the intermediate 
box size. 



from linearity at about 1 cm in all the images. Hence, for q < ~l, the multifractal scaling 
is valid for a shorter length scale range of tank cross section to 1 cm, i.e. ~ 2*. (c) There 
seems to be a further lower cut off of approximately 5 mm, below which the slope of 
J Mij InPij Vs Ib/L for q < 1 goes to zero (See Figure 3(b) at ln(/b/L) ~ —2.75). 
This is expected to be because the box sizes have reached the same order as the plume 
spacings. The geometric structure of the planform does not exist below this length scale. 
To understand the above behaviour, we studied the probability distribution function of 
the measures at various box sizes. The distribution function of the measures is described 
by its moments Zq, and hence is related to r(g) through Zq ^ (Ib/Ly'^'^K The distribu- 
tions of Pij in the standardised form Z — {Pij — Pij) lop^- for Figure |2(h)| where 
indicates the mean and dp. the standard deviation, for three box sizes are shown in 
Figure 01 The measure is distributed normally for the larger box sizes. With decreasing 
box sizes, the low measure and the high measure tails of the distribution deviates from 
the Gaussian. This is clear from the inset in Figure^] which shows the normal probability 
plot of Pij for the intermediate box size; the deviation from linearity shows the difference 
from the Gaussian. All the images in Figure [21 showed similar form and behaviour of 
the probability distribution of the measures. The deviation from linearity for negative 
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Fig ure 5. The / — a sp ectrum for all t he images. .... Figure |2(a)| Figure p(b)| o Fi g- 
ure |2(c)| — Figure |2(d)| x Fi gure |2(e)| Figure |2(f]] A Figure |2(g)| Figure ]2(h)| ★ 

iMeuBveau & SreenjvasanTlllQQr 



Da 



-2.0 -1.5 -1.0 -0.5 0.5 
2.13 2.09 2.05 2.02 2.0 1.98 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
1.96 1.95 1.93 1.92 1.91 1.897 1.886 1.875 1.86 



Table 1. Variation of the Renyi dimension Dq with the moment q for Figure |2(a 



moments in Figure 01 at smaller box sizes seems to be consistent with the deviation of 
the low measure tails of the distribution function in Figure 01 



4.2. The multifractal spectrum 

The f — a curves of all the images in Figure [3 are shown in Figure The error bars in 
the figure are the errors in the slope of the linear curve fits. The value of / at g = is 
the fractal dimension of the 'support' of the measure. As the underlying area over which 
the plumes are formed is a non-fractal surface, /(a(0)) = 2. Figure 13 shows that this 
is satisfied, suggesting a good level of confidence in the calculation. We also note that 
the fractal dimension corresponding to the (information) entropy S — J2f3=oP0^'^EbP0i 
of the underlying process which generates these structures is 1.96. The corresponding 
values Dq for Figure 2(a) are shown in Tabled 

The results of our analysis of the near wall plume structure in turbulent free-convection 
at high Ra are strongly suggestive of its multifractal nature. The analysed images cover 
a wide range of conditions viz., Ra range of about a decade, different test section cross 
sectional areas (10x10 cm^ and 15x15 cm^), absence and presence of a flow through the 
membrane, different large scale flow strengths (about 2.5 times), a wide range of flux (over 
two decades) and show a wide variety of structures having single and multiple large scale 
flow cells with aligned and random structures. The plot of Figure ISjshow that, within the 
errors encountered in the current analysis, all these images have the same multifractal 
scaling of the main plume structure for length scales greater than 0.5cm. 
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5. Discussion and Conclusions 

We have shown that the plan forms of near waU plume structure in high Ra turbulent 
free-convection under varying parameter values have the same multifractal spectrum of 
singularities. Even though the plan forms appear substantially different in their geomet- 
ric form, they have the same non trivial sp atial scaling. The p robabilistic interpretation 
of multifractal formalism, as described bv iMandelbrotI l)l989l) . shows that / and a in- 
directly describe the underlying generating mechanism of these structures in terms of 
the multipliers M for each stage of the process. These exponents are independent of the 
length scales and are unique functions of the multiplier distribution which decides how 
the measure is distributed at each stage. Figure El shows the probability distribution of 
the multipliers in the standardised form, x = {M — M) /cta/- The multiplier values are 
calculated as the ratio of the measures in each box to the measure in its parent box, when 
each box divides into four boxes. The calculations are done for a series of box sizes as 
shown in the legend of the Figure |H| The figure shows that the image has the same form 
of distribution function of multipliers at different length scales, which has an exponential 
tail of higher multipliers. Note that the length scale spans the entire range over which 
/(a) has been obtained. Similar results were obtained for all the images in Figure |21 In 
physical terms, this means that the present images have the same underlying generating 
mechanism. The common generating mechanism that can be identified here is that the 
thin boundary layers of lighter fluid above the membrane becomes unstable resulting in 
the generation of sheet plumes. Interaction of the neighboring plumes due to entrainment 
leads to their merger. New plumes are nucleated in the vacant space due to the merger. 
The whole process is also influenced by the external flow field created by the mean wind, 
whose effect is to align the plume lines in the direction of the mean shear. 

An identical probability distribution of multipliers at different stages of the process 
would show that the multipliers at different stages are not correlated and the under- 
lying process to be statistically self simila r ijMandelbrotl Il989l: IChabra fc SreenivasanI 
Il992l iR-ederiksen. Dahm fc Dowlinglll997l) . The inset in Figure El shows the probability 
distribution function of the multipliers, shown only at two box sizes for clarity. The dis- 
tribution functions coincide at larger multiplier values for the whole multifractal range, 
but are different at low multiplier values. We expect this to be a due to the correlated 
nature of the ends of the structure (low multipliers), which move towards the main plume 
lines due to entrainment. 

Recall that for multinomial processes, f{a) is calculated as the maximum from a range 
of values for any a for e ach moment q. Coii necting this with the thermodynamic analogy 
noted in the literature ijChabra et a/Jll989|) . the common f — a curves might imply that 
the plume structure in turbulent convection is formed so as to maximize the entropy of 
the structure. 

Earlier studies on the mult i fracta l nature of energy dissipation field in turbulent fiows 
by iMeneveau fc SreenivasanI l)l99lh have parallels with the current study. The energy 
dissipation is caused by the velocity gradients at the viscous scales. Plume edges rep- 
resent the major gradients of velocities near the wall in convection. These also are 
the relevant viscous scales near the wall. Hence, it is possible that the near wall en- 
ergy dissipation field follo ws th e plume structure closely. The f — a curves obtained by 
iMeneveau fc Sreenivasa nl lll99lh are shown in comparison with the present curves in Fig- 
ureEl The /(a) curves obtained in our case have a lower spread. Further work is needed to 
clarify the possible connection of these energy dissipation studies to the present analysis. 
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Figure 6. The probability distribution function of the standardised form of the multiphers at 
different box sizes. The inset shows the distribution function of multiphers at two box sizes. 



0.6 




1.8 2 a 2.2 2.4 



0.35 0.4 0.45 I 0.5 0.55 0.6 0.65 
th 

Figure 7. The variation of oidiff = ctmax — ctmin with the threshold Ith for Figure |2(d)| The 
inset shows the f — a curve for the two extreme thresholds 



Appendix A. Effect of thresholding 

The robustness of the calculated multifractal exponents depends on how representative 
the binary image is of the RGB LIF image. Figure [7| shows the variation of adiff = 
Oimax — ctmin wlth the threshold (0 ^ Ith ^ 1) for the image in Figure 2(d) The inset in 
Figure[3shows the f—a curves at the two extreme thresholds. Ith ^ 0.55 produced binary 
images clearly missing in finer details, while Ith ^ 0.45 created images that were noisier, 
with larger plume thickness. Therefore, for this image, the threshold was chosen to be 
0.45. The variation of a^iff for 0.45 ^ Ith ^ 0.55 is about 0.15, of the same order as the 
error in a. Therefore, the change in the a due to the error from the optimum threshold 
is within the error involved in the calculation of a. Further, the inset of Figure [7| shows 
that the threshold does not affect the form of the f — a curve. Hence in all the images, 
the threshold was chosen by visual judgment along with sensitivity analysis. The above 
analysis was for one of the low quality images; for better quality images as in Figure|2(f)| 
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there was negligible variation with threshold. 

We thank G.Vishwanath and Joby Joseph for assistance related to image processing 
and K.R.Sreenivasan for his comments on the draft version. 
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